The role of correlations in the high-pressure phase of FeSe 
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We present a systematic study of the high-pressure FeSe phase performed by means of the first- 
principle electronic structure calculations. Basing on available experimental information about the 
unit cell geometry we calculate the band structure and characterize the related properties during 
their pressure driven evolution. The electronic structure including the hybrid functional B3LYP or 
the Hubbard parameter U for the iron d states lead to the correct semiconducting ground state for 
the hexagonal stoichiometric FeSe within the broad pressure range (up to 30GPa). 



The intriguing class of materials with respect to high- 
temperature superconductivity field is presented by iron 
arsenide based compounds [1-6] which exhibit the crit- 
ical temperatures up to 55 K [7]. The study of these 
materials is complicated by the fact that their magnetic 
and superconducting states are competing at very sim- 
ilar conditions. For this reason the experiments under 
pressure [8] gain special importance since the magnetism 
is typically suppressed by decreasing the volume. In ad- 
dition, the pressure is a cleanest tuning parameter that 
allows to study the electronic structure by ab-initio meth- 
ods. 

Probably the simplest material related to the iron ar- 
senide family is the tetragonal {P4/nmm^ group No. 129) 
/3-phase of FeSe formed by layers of edge-sharing tetra- 
hedrons. Upon cooling the tetragonal phase undergoes 
a slight orthorhombic distortion by lowering its sym- 
metry down to Cmma (Group No. 67) [9, 10]. This 
phase transition takes place within a broad temperature 
range centered at about 100 K depending on a crystal 
size and stoichiometry [9-11]. At rather ordinary tem- 
perature (~ 8K) it becomes superconducting [12], how- 
ever, with pressure the Tc raises amazingly high (up to 
37 K at about 7-9 GPa) [8, 13, 14]. This important infor- 
mation may indicate the direction to search for the new 
superconductors with even higher Tc, which is the main 
focus of research. At first glance, the pressure depen- 
dence of the Tc in FeSe is reminiscent of the supercon- 
ducting dome observed in many unconventional super- 
conductors, such as cuprates, heavy fermions and pnic- 
tides. However in contrast to these systems the vanish- 
ing of superconductivity in FeSe at very high pressure 
is related to a first-order structural phase transition to 
a hexagonal {PG^mmc^ NiAs-type) more densely packed 
phase [8, 15, 16] or its very similar low-temperature or- 
thorhombic modification {Pbnm, MnP-type) [8, 16]. The 
corresponding structures are shown in Figure 1. 

The small amount of NiAs-phase appears already at 
ambient pressures [15, 17] and gradually substitutes the 
superconducting Cmma phase until it fully converts into 
NiAs-type at about 30 GPa. This indicates that at least 



Tetragonal 



NiAs-type 



MnP-type 




Se 



P4/nmm : 129 



Pdsmmc : 194 



* chadov@uni-mainz.de 



FIG. 1. Pressure driven evolution of FeSe phases: su- 
perconducting tetragonal (Cmma) phase is transformed to 
high-symmetry NiAs-type (PGsmmc) structure, which in 
turn transforms into lower-symmetry tetragonal MnP-type 
(Pbnm) phase by little distortion. 



the pressure range of superconducting dome could be ex- 
tended if this structural phase transition is avoided. The 
central aspect in this direction is the knowledge of elec- 
tronic structure of the NiAs-type phase and its related 
properties. Unfortunately there are no systematic stud- 
ies on this system. 

The crucial aspect for FeSe is the stoichiometry. In- 
deed, if the tetragonal low-temperature Fei.oiSe is su- 
perconducting, the small increase of Fe amount (e. g. up 
to Fei.osSe) already leads to a strong magnetic fluctua- 
tions which destroy the superconductivity [10, 18]. The 
analogous situation occurs by substituting Fe with small 
amount of Cu [19, 20]. Many early experiments on non- 
stoichiometric hexagonal phases (e. g. FeySeg) performed 
at the ambient pressure report an antiferromagnetic or- 
der [15, 17]. For thin films the ferromagnetic order was 
reported as well [21]. On the other hand, the stoichio- 
metric hexagonal phase is unstable at ambient conditions 
however even for the lowest pressures the magnetism re- 
veals only in a form of dynamical fluctuations, but no net 
magnetic state was observed. 

The first-principal description of the electronic struc- 
ture under pressure is a difficult task since it requires 
the adequate "total energy-pressure" mapping, or the 
so-called equation of state, based on a full structural op- 
timization at each pressure point. This task is rather 
uncertain even for the high-symmetric structures with a 
single degree of freedom. Indeed, there are different un- 
controlled sources of errors which become especially cru- 
cial for a certain pressure regime ranging from fundamen- 
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tal exchange-correlation inadequacies to the insufficient 
non-sphericity of the one-particle potential. 

Present calculations are based on the so-called CRYS- 
TAL06 code [22] which utilizes the local Gaussian basis 
well-suited for the description of the localized electrons in 
molecules. This formalism can be also applied to describe 
the solid state systems, in particular locally-correlated 
insulators and semiconductors. As a good choice for the 
exchange-correlation potential within this formalism ap- 
pears the so-called B3LYP functional [23] represented as 
the mixture of LDA [24], GGA [25] and the Hartree-Fock 
exact exchange. Since the Hartree-Fock method system- 
atically overestimates the band gap and LDA symmetri- 
cally underestimates it, the mixing coefficients are found 
empirically in order to make use of error cancellation and 
to improve the approach in average for the wide range of 
systems. Since the B3LYP approach typically fails to 
describe itinerant magnetism, our calculations refer to 
the non-magnetic case. Fortunately, this corresponds to 
the experimental evidence which reports the absence of 
magnetic order for the stoichiometric FeSe. 

In the following we perform the full optimization of 
geometry, i. e. without any constraints for the lattice 
parameters and internal coordinates. For the starting 
values we use the experimental structural data [14, 26]. 
Figure 2 (b) demonstrates a very close agreement of the 
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FIG. 2. (color online) (a) Difference of the total energies 
(per atom) of NiAs- and MnP-type structures calculated as a 
function of pressure, (b) The pressure driven evolution of the 
corresponding lattice parameters. The results of optimization 
are displayed by squares. Red circles mark the experimental 
data [14, 26]. 



optimized structures for both MnP- and NiAs-type struc- 
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FIG. 3. (color online) The total DOS curves calculated for 
the optimized NiAs-type (gray shaded area) and MnP-type 
(blue solid line) structures at low (a) and high (b) pressures. 



tures with experiment. As it follows from comparison of 
total energies (Figure 2 (a)) the high-symmetric NiAs- 
type structure is more stable at low pressures, while at 
high pressures the symmetry is reduced to more general 
MnP case. Since this change is caused by a slight dis- 
tortion (see Figure 2 (b)) it does not influence the elec- 
tronic structure as indicated by comparison of the DOS 
curves on Figure 3. For this reason in the following we 
restrict the consideration to more symmetric NiAs-type 
structure. 

As it follows from Figure 3, NiAs-type (as well as MnP- 
type) phase exhibits a semiconducting band gap. This 
finding supports the experimental reports [8]. The cal- 
culated band gap width (about 2.5 eV at ambient pres- 
sure) noticeably exceeds typical experimental values (e. g. 
~ 1.2 eV for the thin films of hexagonal FeSe [27]), how- 
ever such overestimation is typical of all Hartree-Fock 
based methods. 

By increasing the pressure, the band gap gradually 
shrinks indicating the possibility of insulator-metal tran- 
sition at very high pressures (above 30 GPa). Such be- 
havior suggests the strongly-correlated origin of the band 
gap analogical to the situation encountered e. g. in tran- 
sition metal oxides including the known high-Tc super- 
conductors [28-32] which exhibit the pressure driven 
insulator-metal transition accompanied by competition 
of the localized-itinerant electron contributions. 

It is also instructive to probe the correlation-induced 
origin of the gap by applying the alternative approaches 
which account for the local correlations explicitly, as e. g. 
LDA-hU method [33]. On Figure 4 we compare the 
LDA and LDA+U based total DOS curves calculated by 
LMTO (Local Muffin-Tin Orbitals) method [34] within 
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FIG. 4. (color online) DOS curves calculated at (gray 
shaded area) and 11 GPa (blue curve) for the NiAs-type phase 
using LMTO method. Panels (a) and (b) compare the plain 
LDA and LDA+U (t/eff = 7eV) results, respectively. 



scription (t/eff = 0) is indeed inadequate: the minimum 
of total energy is found at about 7 GPa. By increasing 
/7eff, the total energy minimum shifts towards lower pres- 
sures and finally riches the ambient point at [/eff ^ 7 eV. 
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FIG. 5. (color online) Total energy of the NiAs-type phase cal- 
culated as a function of pressure within the LDA+U method 
for UeE = (plain LDA, dashed line) and Uefi = 7 eV (red 
solid line), compared to the B3LYP based calculation (blue 
solid line). 



the so-called PY-LMTO package [35]. Indeed, as it fol- 
lows the band gap could be explained on the basis of 
static approximation to the local on-site electron corre- 
lation, as provided by LDA+U. At the same time the 
plain LDA description which lacks the proper treatment 
of strong correlation delivers a metallic state in agree- 
ment with earlier calculations [36]. The latter is char- 
acterized by a high peaks of DOS at the Fermi energy, 
indicating the instability of the plain LDA solution. The 
band gap is opened due to the strong Coulomb repulsion 
which splits the correlated Fe d-shell into lower and upper 
Hubbard bands. By decreasing the volume the Hubbard 
bands broad, i. e. delocalize and the system metallizes. 

The central parameter of the theory is /7eff = /7 — J, 
where U and J are the effective (screened) Coulomb di- 
rect and exchange interaction potentials. If the J value 
can be relatively easy calculated from the first princi- 
ples and typically does not exceed 1 eV, the screened U 
parameter makes problems. The bare values of U can 
be very high (about 10 eV), however in metals they are 
substantially scaled down due to the intermediate mo- 
bile electrons. The first-principle techniques to estimate 
the screening of U parameter are too imprecise and com- 
putationally demanding in order to be used as practical 
tool. In present case, since the geometry of the system 
is partially known, we can estimate the adequate /7eff 
from the following fit. As it turns out from the pressure 
dependence of total energy (Figure 5), the plain LDA de- 



This huge value of U is actually not much higher than 
typical values needed to obtain the adequate results in 
similar strongly correlated systems [37, 38]. 

Despite the noticeable difference of the band gap values 
(about 1.4 eV within LDA+U and 2.4 eV within B3LYP), 
certain similarity between the Hartree-Fock and LDA+U 
exchange-correlation functionals lead to a similar behav- 
ior of the system under pressure which indicates the Mott 
insulating origin of the NiAs-type phase. 

To conclude, we emphasize that the hexagonal phase 
of FeSe which substitutes the tetragonal superconduct- 
ing phase at higher pressures is characterized as locally 
correlated. These local correlations lead to an insulat- 
ing state which can be classified as Mott insulator. The 
band gap reduces due to increasing electron delocaliza- 
tion with pressure up to the insulator-metal transition 
which occurs roughly above 30 GPa. These results rea- 
sonably correspond to the experimental studies of FeSe 
and related systems under high pressure. 
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